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0.  FOREWORD 

The  work  done  during  the  period  covered  by  the  present  FINAL  Scientific 
reports  pertains  to  the  following  topics  related  to  the  application  of 
functional  analysis  in  fluidmechanics: 

1)  Variational  formulation  of  non-linear  boundary  value  problems; 

2)  Closed  splines. 

The  present  Report  is  divided  in  two  parts. 

The  first  part  describes  more  or  less  concisely,  according  to  the  needs, 
the  work  which  has  already  been  extensively  reported  in  available  publications 
It  consists  of  three  chapters  dealing  with  : a)  A class  of  non  linear 
boundary  value  problems;  b)  Applications  of  closed  spline  functions; 
c)  Numerical  sperimentations  with  closed  splines. 

The  description  of  the  nature  of  this  work  and  of  the  publications  it  has 
originated  is  given  in  what  follows. 

As  a natural  development  of  the  research  program  on  applications  of 
fun  ■' ionai  analysis  to  fluid-mechanics,  the  main  objective  of  the 
next  stage  of  the  investigation  was  clearly  to  be  identified  with  the 
extension  of  the  results  found  for  linear  problems  to  as  large  a class  of 
r.on-linear  problems  as  possible. 

A first  positive  step  in  this  direction  has  been  accomplished. 

The  following  prob Lem  has  been  investigated  and  solved:  variational  for- 
mulation of  the  non-linear  problem  A (u)  = f(uj  where  A is  a positive 
definite,  formally  self-adjoint  operator  in  a suitably  defined  Hilbert 
space  H of  the  functions  ( ^ ) and  | is  a Frechlt.  differentiable  non-linear 
operator  defined  in  H. 

A first,  version  of  a paper  describing  this  work  and  co-authored  by  the 
principal  investigator  and  by  Dr.  C.  G0I.1A  lias  been  presented  at  the  Third 
Congress  of  the  Italian  Association  for  Theoretical  and  Applied  Mechanics 
IAFMKTA)  held  in  Cagliari  (Italy)  during  the  month  of  October,  I97e  Ref.  ll 
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In  the  preprint  of  this  paper  acknowledgment  to  AFOSR  support 
vsas  inadvertedly  omitted  and  only  the  sponsoring  received  by  the  second 
author  was  acknowledged.  This  error  will  be  corrected  in  the  second 
version  of  the  paper  (see  later). 

Further  studies  on  the  subject  have  led  to  a much  improved  formulation 
and  solution  of  said  problem.  This  has  resulted  in  the  paper: 

1.  L.G.  NAPOLITANO,  C.  GOLIA : " Dual  Variational  Formulation 
of  a non-Linear  Boundary  Value  Problem" 

which  will  be  submitted  for  publication  in  an  International  Journal. 

Since  the  first  version  of  the  papar  is  available  [»'■  i]  and  the 
final  version  presents  only  formal  improvements  of  the  first  one,  only 
a concise  description  of  the  results  needs  to  be  given.  T-is  will  be 
found  in  Part  I of  the  present  Report,  chapter  one. 

The  activity  concerning  closed  splines  has  been  concerned  with: 

i)  working  out  applications  of  the  closed  spline  functions  introduced 
by  the  principal  investigator  jjlef.  2~\ 

ii)  carrying  out  numerical  sperimentation  with  these  closed  splines 

iii )  defining  and  studying  the  properties  of  a new  class  of  closed  splines: 
the  Hermite  closed  splines 

The  definition,  characterization  and  properties  of  closed  splines  were 
described  in  the  Final  Scientific  Report  of  the  previous  Grant  AFOSR 
76-2889  [Ref.  2 ] . These  new  classes  of  splines  yielded  interpolating 
functions  for  closed  curves  and  possess  all  the  properties  of  spline- 
functions.  Conventional  approaches  still  largely  use  conventional 
splines  to  interpolate  closed  curves:  the  analysis  reported  iri  [ 3 ] illus- 
trates the  many  shortcomings  of  such  an  approach  and  shows  how  markedly 
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superior  is  the  use  of  closed  splines.  To  clearly  substantiate  these 
theoretical  predictions  it  was  felt  appropriate  to  work  out  some  examples 
of  applications. 

The  problem  of  interpolating  significant  airfoil  shapes  was  considered 
and  results  obtained  from  classical  and  closed-splines  formulations  were 
compared. 

These  examples  of  applications  of  closed-splines  are  described  in 
Part  I of  the  present  report,  chapter  two.  They  were  included  in  the 
paper : 

2.  L.G.  NAPOLITANO,  V,  LOSITO  " The  closed  splines  functions" 

which  has  been  accepted  for  publication  in  the  "Journal  of  Computational 
Physics" . 

Once  the  theory  of  closed  spline  functions  had  been  developed,  the  next 
stage  of  analysis  was  addressed.  It  concerned  a number  of  questions 
related  to  the  optimal  use  of  closed  splines  in  a computational  scheme 
and  called  for  extensive  numerical  sperimentation . 

This  phase  of  the  work  is  completed  and  its  results  are  inscribed  in  the 
following  paper  which  has  been  accepted  for  presentation  at  the  Congress 
of  the  Italian  Association  for  Aeronautics  and  Astronautics  (AIDAA)  which 
will  be  held  in  Milan,  during  the  last  week  of  Septebmer  1977: 

3.  L.G.  NAPOLITANO,  V.  LOSITO  and  A.  VITIELLO  "Numerical 
sperimentation  with  Closed  Splines  " 

For  this  reason  only  a concise  summary  of  theprobleras  addressed  and 
of  the  results  obtained  will  be  given  in  Part  I of  this  report.  Chapter 
three. 


The  second  part  of  this  report  describ*'  n extenso  the  third  development 
of  the  research  work  on  closed  spline,  i.c.  the  one  related  to  Hermite 
closed  splines. 
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PART  I SUMMARY  OF  WORK  ALREADY  PUBLISHED 

1.1.  INTRODUCTION 

As  mentioned  in  the  Foreword,  this  first  part  of  the  report  describes 
concisely  work  which  is  already  available  or  will  shortly  become  available 
in  published  form. 


1.2.  DUAL  VARIATIONAL  FORMULATION  OF  A CLASS  OF  NON  LINEAR  BOUNDARY 
VALUE  PROBLEMS. 


The  following  class  of  non-linear  boundary  value  problems  has 

been  considered,  where  A is  a positive  definite,  formally  self-adjoint 
operator  in  a suitably  defined  Hilbert  space  H of  the  functions  ( u ) and 
f is  a Frechet  differentiable  non-linear  operator  defined  on  H. 

First  the  conditions  on  f which  allow  a strong  variational  formulation 
of  the  problem  have  been  determined.  It  has  been  found  that  these  conditions 
reduce  to  the  symmetry  of  the  non-linear  mapping  f (u). 


Subsequently  the  hybrid  strong  variational  formulation  has  been  derived  from 
which  the  dual  variational  formulations  have  been  obtained  as  particular 
cases  £ see  Refs  (j>  ) , c«S  originated  under  a previous  grant  for  the  meaning 
which  is  ascribed  to  hybrid  functionals  in  the  subject  context/. 

Finally,  upper  and  lower  bounds  for  the  exact  solution  and  error  estimates  for 
the  approximate  solutions  have  been  found. 


The  approach  used  is  based  on  the  search  of  an  hybrid  functional  by  means 
of  integration,  in  a cartesian  product  of  two  suitable  Hilbert  spaces,  of 
a gradient.  Imposition  of  the  integrability  conditions  and  of  appropriate 
stationarily  conditions  leads  to  the  requirements  on  A and  f that  must  be 
satisfied  for  a strong  variational  formulation  to  be  possible.  The  subsequent 
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A number  of  numerical  techniques  for  the  solution  of  planar  inviscid  flow 
fields  around  airfoils  require  knowledge  of  the  parametric  equations 

'X.(yi')  j of  the  airfoil's  profile  in  terms  of  its  arc  length  2T. 

These  equations  must  be  found  by  interpolating  the  data  J representing 

the  ordinates,  of  the  upper  and  lower  profile  at  a number  of  stations  x.;  . 

The  arc  length  2T  is,  of  course,  not  known  a priori  and  the  problem  has  to 
be  solved  iteratively.  A possible  procedure  is  as  follows. 

Assume  as  curve  parameter  the  length  w measured  along  the  polygonal  con- 
necting the  assigned  profile  points  to  obtain,  from  the  data  z^  = z (x.) 
two  sets  of  values : 

where  w^  = w(P^) 

From  these  two  sets  construct  the  two  interpolating  functions  x(w) ; z(w)  and 


compute 


.y  V , ) f fe[z\z+  /gIk 

t J l (#Ji W \e(i4/j 


2 7^ 


The  quantity  w^^  £ which  would  give  the  profile's  arc  length  if  x(w), 
z(w)  where  the  exact  parametric  equations  of  the  profile  in  terms  of  w J 
is  taken  as  the  new  curve  parameter  and  the  cycle  is  repeated.  From  the 
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new  sets  of  values 


*/ 77/' 


two  new 

constructed,  a new  curve  parameter 


interpolating  functions  k /Vl''7;  zfv"]*- 


are 


U/y 


determined  and  so  on. 


<*■>-  ^ 


The  iteration  terminates  when  / w^'-  '7  m-4LX  is  smaller  than  a pre- 


(r) 


assigned  small  quantity  for  then  the  r-th  iterate  w pratically  coincides 
with  the  arc  length  ( Z ) and  the  functions  X £ IV  J ' Z £ ikewise 
coincide  with  the  interpolating  functions  X (Z)  ^ Z(jc-  ) • 


At  each  cycle  the  prpcedure  requires  the  construction  of  two  interpolating 
functions  x[vu’l  z[tvLiT]  . In  the  examples  to  be  discussed  here 
conventional  and  closed  cubic  splines  (q=2)  with  several  values  for  n are 
considered. 

The  interpolating  functions  based  on  conventional  splines  are  denoted 
without  any  subscript  and  those  based  on  closed  splines  with  a subscript  (c). 
To  simplify  the  notations,  primes  will  always  denote  derivatives  with  respect 
to  the  pertinent  parameter. 


According  to  the  results  of  the  theory,  £ see  Refs  (_ 2 ) ,(3)]for  any  (n)  and 

at  each  cycle:  x'  , x,(,  z',  z ' ' will  be  continuous  everywhere,  x'  and  z' 
o c c c 

will  be  discontinuous  at  the  closure  point,  x’’  and  z’’  will  be  continuous 
everywhere  but  will  vanish  at  the  closure  point.  Hence,  at  this  point, 
the  conventional  splines  will  yield  an  infinite  radius  of  curvature  R. 


To  facilitate  comparison  with  exact  values  a family  of  elliptic  profiles 
is  considered,  which  is  defined  by  the  equations  x = 5(l-cos  9) , z = jV2  sin  9 
where  9 is  the  polar  angle  and  0.5  £ ? - 1 is  the  thickness  ratio. 

The  closure  point  is  the  point  corresponding  to  9 = 0 and  the  other  points 

correspond  to  constant  increments  in  0 • 


The  question  of  the  optimal  location  of  the  points  is  an  interesting  one 
and  becomes  very  important  for  low  values  of  ^>  . It  will  be  addressed  in 
the  next  paragraph. 


8 


The  discussion  will  be  limited  to  the  final  interpolating  parametric 
curves  x ( ) , z ( 2"  ) as  in  all  cases  considered  the  iteration  converged 

very  rapidly  (at  most  two  cycles  were  needed). 

Due  to  the  symmetries  of  the  airfoils  and  to  the  choice  of  the  points  , the 
discontinuity  in  the  drivative  z'  at  the  closure  point  turns  out  to  be  zero 
whereas  that  of  the  derivative  x'  is  equal  to  twice  the  value  of  x'  at  Q = 0 . 

The  values  of  the  discontinuities  x'  at  the  closure  point  for  different 
values  of  (n)  and  for  elliptic  airfoils  (^>=  1,  .75,  .5)  are  shown  in 
Tab.  1.  Obviously  they  do  decrease  as  n increases,  but  the  point  here 
is  that  a much  better  result  (vanishing  discontinuities)  can  be  obtained  at 
no  cost  (i.e.  without  having  to  increase  the  number  of  subintervals)  with  the 
use  of  closed  spline  functions.  Tab.  I also  shows  that  things  get  worse 
as  ^ decreases:  the  same  trend  will  be  found  in  all  other  features. 

Values  of  the  derivatives  dx/d  f and  dz/dT  are  shown  in  Tab.  Ila  through 
lie  for  the  same  three  values  of  ^ and  n=8.  Due  to  the  symmetry  of  the 
profiles  only  values  for  0 — 0 £ are  listed:  the  first  column  gives 

the  exact  values,  the  other  two  columns  give  the  values  obtained  with  the 
closed  spline  and  the  classical  spline,  respectively. 


3 anticipated  in  Refs  (^2)  , ^3)  : 

i)  the  dis conti nuity  at  the  closure  point  is  associated  with  poor 

accuracy  in  its  neighborhood.  Thus,  since  wv  0 whereas  Sz'  =0  at 
the  closure  point,  the  accuracy  in  x'  is  worse  than  in  z*  in  its 
neighborhood.  Tables  II  show  that  both  accuracies  decrease  as 
decreases  and  that  the  results  yielded  by  the  closed  splines  are  far 
superior  to  those  given  by  the  conventional  ones; 

ii)  the  discontinuity  at  the  closure  point  influences  the  behaviour  of 

the  interpolating  function  throughout.  Thus,  for  instance,  the  inter- 
polating functions  based  on  classical  spline  fail  to  exhibit  the  symmetry 
properties  dx/d  Z («/^  +»)  * dx/d  Z('f  i^~  & ) and  dz/dT  ( /T/2  + Q ) = 

- dz/dT  (T/o  - B ) , which  are  instead  satisfied  by  the  interpolating 
functions  based  on  closed  spline; 
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iii)  as  one  moves  away  from  the  closure  point  the  "smoothness  properties" 
of  the  interpolating  functions  based  on  classical  spline  tend  to 
approach,  at  a rate  depending  on  the  curve  to  be  interpolated,  those 
of  the  interpolating  functions  based  on  closed  splines. 

Exact  and  approximate  values  of  the  radius  of  curvature  R were  also  analyzed 
for  ^ = 1 and  0.5  with  n=8.  The  analysis  confirmed  that: 

l)  the  classical  spline  interpolating  function  is  vastly  off  in  the 
neighborhood  of  the  closure  point,  whereas  the  closed  spline  gives  there 
very  satisfactory  values; 

<-)  the  performance  of  the  classical  spline  approaches  that  of  the  closed 
spline  sufficiently  far  away  from  the  closure  point. 


10 


TAB.  [ 

Discontinuities  ax'  at  the  closure  point  for  different  values  of 
n when  elliptic  airfoils  are  approximated  by  means  of  classical 
spl i nes . 


.2322  .1139 

.3306  . 1646 


.5546 


.2444 


1 . U . NUMERICAL  SPERIMENTATION  WITH  CLOSED  SPLINES. 


As  discussea.  in  the  previous  paragraph  the  performance  of  closed-splines 
interpolating  functions  tended  to  become  worse  and  worse  as  the  percent ua] 
thickness  of  the  elliptical  airfoils  decreased. 

This  tendency  is  to  be  ascribed  to  the  fact  that,  as  ^ decreases,  the 
curvatures  at  the  leading  (and  trailing)  edges  of  the  airfoil  increase 
whereas  that  of  the  central  part  decreases.  Hence  there  was  a clear 
indication  that  an  uniform  polar  angle  ' spacing  of  the  data  points  is 
not  always  adequate . 


The  optimal  distribution  of  a prescribed  total  numer  (n)  of  points  P 
constituted  the  first  objective  of  numerical  sperimentation  with  closed 
spline . 


The  obvious  speculation  was  that  the  points  P^  had  ti  be  "concentrated"  in  the 
large  curvature  region.  Numerical  sperimentation  confirmed  this  speculation 
and  showed  a marked  improvement  in  the  accuracy  of  the  interpolating  functions 
with  the  use  of  cubic  closed  splines. 

These  encouraging  results  point  to  a very  promising  area  of  further  theoretical 
research,  namely  that  of  optimal  location  of  data  points.  Some  preliminary 
suggestions  on  possible  theoretical  approaches  are  given  in  the  paper  mentioned 


eword . ^Ref  (.Wj. 


in  t.;  e Fore' 

CC 

The  relation  remedy  is  naturally  applicable  only  when  one  can  prescribe 
arbitrarily  the  position  of  the  points  P. . 

Whet,  this  is  not  the  case,  i.e.  when  both  the  total  number  of  points  and 
their  locations  are  given,  the  alternative  may  seem  to  be  an  increase  in  the 
degree  (?)  of  the  closed  spline. 


Extensive  numerical  sperimentation  proved  however  that^in  general,  this  is 
not  so.  It  turned  out,  indeed, that  for  each  number  ( n ) of  points 
P.  (and  given  distribution  of  them)  there  is  an  optimal  value  qQ  of  the 
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degree  of  the  spline  that  gives  the  ''best"  approximation.  Or,  equivalently, 
for  any  given  degree  of  the  spline  there  is,  for  a given  distribution ^ and 
optimal  value  n^  of  point  which  yields  the  "best"  approximation. 

The  optimal  value  yi0 ^ ifj  ) increases  with  q.  Thus,  els  stated,  it  is  a 
fallacious  illusion  the  one  that,  for  a given  fixed  number  of  points 
and  a given  distribution  of  them,  hopes  to  increase  the  accuracy  of  the 
interpolation  by  increasing  the  degree  of  the  spline.  This  will  be  the 
case  if  and  only  if  the  initial  degree  of  the  spline  is  less  than  the 
optimal  one  corresponding  to  the  given  number  (n)  of  points. 

This  being  the  state  of  affairs  the  true  problem  is  then  how  to  improve 
on  the  accuracy  of  the  approximation  when  both  the  number  of  points  (nad 
their  distribution)  and  the  degree  of  the  spline  (at  its  optimal  value) 
are  fixed. 

Also  for  this  problem  some  preliminary  suggestions  are  offered  in  the 
paper  mentioned  in  the  Foreword  ^ Ref.  ( 
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PART  II 


HERMITS  CLOSED  SPLINE  FUNCTIONS 


*>y 


L.  G.  NAPOLITANO 


1.  INTRODUCTION 

In  computational  aerodynamics  it  is  often  necessary  to  solve  interpo- 
lating problems  related  to  "airfoils",  i.e.  to  closed  curves. 

Using  classical  spline  functions  to  interpolate  airfoil's  ordinates 
prescribed  on  a finite  set  of  point  is  unsatisfactory,  and  on  many 
important  accounts  Cl]  and  a markedly  better  approach  is  afforded  by 
the  closed  spline  functions  li]. 

A similar  situation  arises  when  the  data  to  be  interpolated  represent 
values  of  the  functions  and  of  its  first  derivative  at  a given  set  of 
points . 
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Consider  indeed  the  problem  solved  by  classical  Hermite  spline 
functions . 

Given  a closed  interval,  reduced  through  suitable  normalization  to 
0,lj^  (n+l)  points  fr  ) with : 

o = «■<  « ^ * %n.  < i 

and  (2n)  real  numbers  (ri);  (ri')  ( / < L < Yl+!)  the  interpolating 

Hermite  spline  function  of  order  Cj  < YL  i / corresponding  to  the 

(n  + l)  points  T?*  and  to  the  sets  \ 1 ? is  the  unique 

-V  /ll-  -y  L J J L L J 

function  £ f-f  /,  Oj  I J which  solves  the  following 

minimization  problem: 

f ynCn  f [ f^\z) iVr 

j 

[ J '(*)--  K' 


where 


denotes  the  first  derivative  of 


p) 


and 


U\o,l] 


vs  the  Hilbert  space  of  real  functions  ? ( ) defined  on  f 0 I 7 and 

i ^ y oL*i)  J 

having  a square-integrable  q-th  derivative  K " J • 

The  space  r[o,i  7 quantifies  the  nation  of  "degree  of  smoothness" 
of  the  interpolating  curve. 

The  Hermite  interpolating  splines  belong  to  the  2.  + 0 - dimen- 

ior-i  cubspace  S H {.&,!]  °f  real  functions  i 

defined  on  L^/l] and  such  that  f*J 


*V 

a)  A(jcji-S  a poiinomial  of  degree  (2q~i)  in  eacli  of  the  open 
intervals  ] J (l^L$  Vi ) 

<-v  i J arid  i tv  iirsi  (2q-3)  derivatives  are  cwntin^Oua 


v.  j 
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n 


i 


c)  the  derivatives  of  d (£)  f rom  ordersVto  (2q-l)  vanish  at 

the  points  X.  = 0,  'C  - 1 . 

' Iff  I 

Such  an  Hermite  spline  is  given  by 


V 

-i  Cz)~ 


J 


~ c 


as  £±3 

j J!  (iij-iJJ 


'1.2  ‘ 


where : 


tr-Zi J+  -f  d.3) 


(r-r,j+- 


o 


it  (J-  4 0 


and  the  2 (n+l)  coefficients 


i.J 


satisfy  the  q equations: 

k~l 


(l.U) 


ntl  r - ^ , 7l  r \k~‘  7 ^ 

Z [A&)  4 /'  ° 

d~'  [o  < k z ? 

l h.-i 

(with  the  convention  that  R.  X-  = 0 for  k = 0). 

By  imposing  the  2(n  + l)  additional  requirements 

iZr; j Z'teyt;  ; ' 

the  |_q  + 2 (n  + l)  J coefficients  are  uniquely  determined 

2 "]  and  so  it  is  the  resulting  Hermite  interpolating  function  (P) 


(1.5) 


L 
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solving  problem  (l.l) 


If  the  data  to  be  interpolated  correspond  to  a closed  curve  then:  i,1—  ‘<^4/ 

. / . { 

(closure  point);  V - 1 Z.  ~ Z . and  Z"  is  to  be  inter- 

preted  as  the  curvilinear  coordinate  of  a sufficiently  smooth  and 
regular  closed  contour  C normalized  with  respect  to  its  length  jfl  7 • 

If,  to  interpolate  these  data,  one  uses  the  classical  Hermite  function 
& [z  )=  ^ (_  whe  re  S^rjis  defined  by  (J  , 2. J through  (1.5) 

with  %,-t  then  O'f'clwill  indeed  have  the  above  mentioned 

spline  properties  in  each  of  the  (n)  open  sub-intervals  7 ^ ~ M * 

However,  at  the  closure  point  P,  - 

> 1 

i)  the  conditions  c ^ + ^ guarantee  only  the  continuity 

of  the  function  and  of  its  first  derivative;  the  other  derivatives,  up 
to  the  order  (q-l),  will  not  be  continuous  ^unless  the  location  of 
points  "Fj  and  the  data  sets  J satisfy  particular 

s^punetry  conditions  J , their  discontinuities  being  uniquely  determined; 

£"  upon  the  uniqueness  of  the  solution  of  problem  (l.l)J  J 

ii)  the  subsequent  derivatives,  up  to  the  order  (2q-3),  are  continuous 
but  van  1^1  identically  £ upon  the  propertis  (b)  and  (c)  of  3 [Z  J 


iii)  the  discontinuity  of  the  remaining  derivatives  are  likewise 
uniquely  determined:  that  of  the  (2q-2)-th  derivative  being  equal 
to  1L-  o and  that  of  the  (2q-l)-th  derivative  being  equal  to 

< l - L.,  >• 


Thus,  as  for  the  problem  considered  in  [*1  J , the  use  of  conventional 
Hermite  spline  functions  for  closed  curves  presents  a number  of  shortcomings. 
Unwanted  sicontinuites  are  introduced  at  the  closure  point  and  the  vanishing 
of  the  derivatives  from  orders  q to  (2q  - 3)  may  be  too  penalizing, 
especially  for  low  values  of  q.  In  addition,  always  because  of  point  i) 


21 


above,  the  -terpolating  function  does  not  belong  to  H1^  (C).  Hence, 
its  "degree  of  smoothness"  is  not  characterized  and,  perhaps  more 
critically,  <5"  does  not  satisfy  any  minimum  problem. 

These  shortcomings  are  eliminated  if  one  uses  the  notion  of  closed 
splines  introduced  by  the  author  in  £"lj  and  constructs  a new  class  of 
such  splines,  the  Hermite  closed  splines,  constituting,  as  the  others, 

Q,  ' ) • fl  ' V ••  •'  - •'  'T  1.  . o jj  1 X.  O ^4  0.  O 04  V>  VW  uOun  • 


The  approach  to  be  used  to  construct  and  study  Hermite  closed  splines 
is  based  on  the  abstract-space  spline  theory  detailed  ’ n ^ 

perhaps  appropriate  tojustify  its  use. 

One  might  argue  that  all  one  needs  is  the  introduction  of  (q-2)  additional 
parameters  to  impose  the  missing  continuity  of  the  derivatives  from 
orders  2 to  ( q - l).  This  speculation  is  substantially  correct.  The 
resulting  interpolating  spline  would  have  a total  of  [ ^ + 
parameters  and  would  essentially  coincide  with  what  we  call  a Hermite 
•jiuseu  interpolating  splir.i.  However  one  would,  reel  at  loss  as  to  how 
actually  to  proceed,  and,  more  fundamentally,  questions  related  to 

exi  4 en-""  • d ■ ■ - w i *>■•••»  relating  functions  would  be  left 

unanswered.  Furthermore  it  would  be  difficult  to  ascertain  and  eventually 
establish  *he  eauivahnet  . •>•...  i • poialion  problem  and  a 

minimization  prop  Lem  to  identify  the  subspace  generated  by  these  "splines"^ 
to  describe  their  extremal  properties  and  so  on. 

*'.r  ? u tiora-  j by  contrast.,  are  most  simply  resolved  by  using 
abstract  space  snline  theory. 

,-.o  iii  j < j .to  facilitate  a iuore  -p  <*i-  c mp  (‘pension , the  definition 

characterization  and  properties  of  Hermite  closed-splie  functions  are 
■first,  stall'd  without  proof  in  paragraph  (2). 

Existence  uniqueness  characterization  extremal  and  other  relevant 
properties  are  proven  in  paragraph  (d),  wherein  for  completeness'  sake,  the 
pc  i u4iicnt  iieeuen  result*  f”T  th*  '''Vhc*'  -.a  or  *h*?ry  of  spline-functions 


k 
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are  also  concisely  recalled. 

As  this  paper  employs,  with  few  variants,  the  same  approach  developed 
in  [ 1 ] it  will  be  rather  more  concise  and  reference  is  to  be  made 
to  [l  } for  missing  details. 

2.  DEFINITIONS,  CHARACTERIZATION  AND  PROPERTIES  OF  HERMITE  CLOSED  SPLINES 

a)  Defi  nit  ions 

Let  C be  a sufficiently  smooth  closed  contour.  Denote  by  £"  the  curvi- 
linear coordinate  along  ' measured  free  ur,  art.!  ! rary  initial  point  P^ 
and  normalized  with  respect  to  the  length  of  C,  so  that  O d t d , 

to  be  referred  to  as  the  closure  point  on  the  contour,  is  characterized 
by  either  value  2"  — O ' Z - d . 

Consider  (n)  arbitrary  successive  points  5JfJbn  C:  let  Z-  be  their 

curvilinear  coordinates  with: 

o ~ r,  < r2  < ?3  ... < «u  i 

/ 

and  prescribe  (?n)  r*» ai  numbers  c ‘ V '. 

^ J ir 


The  Hermite  closed  interpolating  spline  <5~of  degree  (q)  corresponding  to 
the  (n)  triplets  ( 2Tj  * j ± <■  - YL.  ) is  defined  as  the  unique 

element  of  H*!  (C)  such  that: 
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In  order  for  <T£  hV)  to  be  such  a Hermite  closed  interpolating  spline 
it  is  necessary  and  sufficient  that: 

a)  6"  be  a polynomial  of  degree  (2q-l)  in  each  open  sub-interval 

1 r-y  


b)  he  continuous,  on  C,  together  with  its  first  (2q  - 3)  deri- 


vatives, i.e. 


Lkl(o*)-  ^ k>  (i ')  ( Y k 6 Lc -3] 

o'%VMJ^/7vy  ; *ieC*.*l  J 

c)  0"  be  such  taht 

<r(t; > K-  ; <?’(*;)=  *i  ) v 1 6 

The  set  of  functions  ^(Z)  & H 0- ) satisfying  the  first  two  conditions 

constitutes  a subspace  .S  [~the  space  of  Hermite  closed  spline  functions 

Lc  f 

corresponding  to  the  set  J Z^- 1 J °T  dimension  (2n). 
b ) Characterization 

I 

An  element  of  Hq(0)  belong--  t;  the  :ul..;  »«•-  O o:  Hermite  closed  splines 

ly  , Ut 

corresponding  tcfc^fiP  it.  is  representable  as 


[z)*k  9 „ £./  LL-zj 

j-°  J J!  ^'t-ux-crr 


2<?-’  . / , 2 “ 2 7 

t + ^ / ’ 

where  the  (2n)  coefficients  La.  and  the  ( ?q ) cccfficicnts  IB  satisfy 

/ i I i 

the  following  (Pq)  equations:  ' 


2.2) 


2h 


2 A 1*0 

i-  I 


J 


176,.,,,  £s¥]j 

:c  Tf^V  U Cc~l)1  (t-rOi  J 


(2.3) 


(2.4) 


The  functions  (Jt'Zi J is  defined  by  equation  (1.3):  it  is  continuous  on  C 
together  with  its  firsf(k  - l)  derivatives  and  its  k-th  derivative  is  dis- 
continuous at  C~  2^  , the  discontinuity  being  equal  to  ( fe !)  [*]- 

The  (2q~l)-th  and  (2q-2)-th  derivatives  of  s ( ) are  given  by: 

♦ B 2-  ♦ 2-  M 

'2Sj-z  ]2Sf-,  o = > 

Hence : 


- nl  tJ-  &- 


i)  the  (2q-l)-th  derivative  of  s ( Zr  ) is  piecewise  constant  in  each 

L~-  j 

open  sub-interval.  Specifically  it  is  equal  to  ( G ■+  7>  X / ) 

T ( ' . \ O A — 1 * 

ln  J £ C;  . , L , U - L - % ~l ) and  to  p in  7 •/  r . Thus  the  (n) 

/ , ' nv  *- 

coeffi  cientc,  ^4*  represent  tiie  values  of  the  discontinuity: 


(*«!  -o 


at  the  Points  . Notice  that  for  C ■=/  it  is  Z,  ^ O y Z,  = / ; thus 

S'  -t^-P  --  ^ 

x 'a^-/  »2.^-/ 

ii ) The  (2q-2)-th  derivative  of  s ( 2T)  is  piecewise  linear  in  each  open 

sub-interval.  Specifically,  it  is  equal  to 

L 


z-n)  + At  ] 

~-l  - '2S-' 

in  ti  /l_  ( I 6 t and,  on  account  of  eqs  (2.3)  and  (2.U) 


ft- • P..7  - 5,  i t- ( 
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for  r = 1,  to 

R -r  G [c-  l ) 

, - T>. 

in  J C ti , I L ' The  discontinuities  of  j,  at  the  points  I L 

are  thus  given  by: 

?r'2,c?o=  A':  c*i£fn-) 

s^’cv,)  = (pa,.,t^  '8,.,'  /t- 

Hence  the  coefficients  represent  the  discontinuities  of  the  (2q-2)-th 
derivatives  of  s ( Z ) at  the  points  Vc  [i*i  ± *?]• 

Notice  how  the  closure  point  ( R ) is  in  no  way  differentiated  from  the 
other  (interior)  points. 

Given  (2n)  arbitrary  real  numbers  £4 ^ fcj  , thre  is  a unique  Hermite  closed 
spline  corresponding  to  the  set  Zf  £ such  that: 


(2.5) 


Conseauentlv , the  system  of  2(n  + q)  equations  (2.3),  (2.4),  12. b). 


' f ) . w w . 


2 <n  + q)  coefficients  appearing  in  the  unique  interpolating  Hermite  spline 
fun oti  which  solves  the  minimum  problem  (2.1). 


c ) Properties 


5 r r y ■■  ••m  nl  4,4.  ^ and  apy  fraction  P £ //(^,)  it  ;s: 

dfVw rfr-t  hn)-XL  frj 

J t:i  ' 


(2.6) 
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In  particular: 


l-  jfe mo,*]] 


The  minimization  problem  (2.1)  solved  by  ^ [zj is  a particular  case  of 
the  following  more  general  extremal  properties. 

Given  (n)  arbitrary  but  fixed  xj  is  the  unique  element 

of  SUc  such  that  Vl7-)='Z-$[Z):t.then, 

f\  ) Fo'fc.  cu^.  I; 

--  ^ fc  £ ^ 


is  the  unique  element 


f[^V-fV]^z 


(2.8) 


and  any  other  element  CTof  vS  having  this  property  differs  from  (S’  by 
a constant; 

J^for  any  ^ 6 -S^c  ■' 

^Telc  =rn^  l ] dt  {!9) 

and  5”  is  the  unique  element  of  T having  this  property. 

If  the  r.h.s.  of  equations  (2.h)  is  denoted  by  , the  discontinuity  Si  OP,) 
of  the  k-th  derivative  of  ^ at  the  closure  point  P,  is  given  by 


S ) - - L 


itf-k-i 


l 
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Hence  of  the  (2q-l)  homogeneous  conditions  given  by  equations  (2.1*)  only 
the  last  (2q-2)  correspond  to  the  required  continuity  conditions  at  the 
closure  point.  The  remaining  one  (for  » l) , which  reads: 


& B - 

2 ^ -f  c - ' 


A;  zc 


A:  » C? 


(2.10) 


1 ■:'T' 1 y expresses , as  equation  ^2.3),  obvious  congruence  conditions  for 
the  discontinuities. 

Indeed,  if  h(z)is  any  function  discontinuous  at  the  point  fj  and  if  we 
let: 

*u  - ^ )-  L fc'J 

then, 

A L*i\,)-h(TC)*  A-Ltlk(V,H.) 

so  that.,  by  summing  for  (i)  from  ( i ) to  (n)  and  considering  that, 

by  the  definition  of  the  closing  point,  P =P  • 2"  - 2 * = 0*  one  obtains 

J J V + r 

tnat  the  discontinuities  b h ( Pj J must  satisfy  the  "congruence"  con- 
dition: 

Z 4 . k * S I (ft  Jr  C? 

< - / b - / . 

which  is  really  checked  to  reduce  to  equation  (2.3)  for  h (r)  = ^ fe ) 

/ £ /t  ( Pj  J - J--  J J {/To  equation  (2.10)  for  h ( 2T ) = ££/)  =/(.£,: 


and 


ince 


it  is  accounted  for  the  fact  that: 


j 


~'J  (Zq  ~/J 

where  3 . is  the  constant  value  of  i 

subinterval  J Tc  ? L 4 , L 


in  the 


“d  * / ; z - O. 
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3.  HERMITE  CLOSED  SPLINES:  PROOFS  OF  EXISTENCE  UNIQUENESS  AND 

CHARACTERIZATION 


Proofs  of  the  statements  made  in  the  previous  paragraphs  hinge  on  the 
Hilbert  space  formulation  of  spline-functions  theory  detailed  in  2 J . 

For  the  sake  both  of  self-completeness  and  of  ready  reference  the  basic 
results^as  directly  specialized  to  the  present  nee ds,  are  briefly  summarized. 
A slightly  more  general  formulation  was  given  in  flj  : greater  details 
are  to  be  found  in  Ul 

a. 

If  C is^Tufficiently  smooth  and  regular  closed  contour,  consider  the 
two  Hilbert  spaces  with  their  standard  inner 

products  [ denoted  by  <f  , ^ and  j > ^ o J and  the 

linear  continuous  operator  (c)~>  (f-J  where  denotes  the  q-th 

derivative  with  respect  to  the  previously  defined  curvilinear  coordinate 


2T  <5  L 0/  / J , normalized  with  respect  to  the  length  of  C 
Finally  let  (2n)  functionals  be  defined  by: 

< K;,  = *(ji )*  xe  X-tfXc) 

< ; x>^  = /i'Ci)-  t’i  i < l ± n. 


(3.1) 


0=  r,  < * 


< r c 4 


c_ 

By  letting  Zc  R j a finite  dimensional  subspace,  and 
Z * , j ,tz  f J ) £ Z. 


equat 


ions  (3.1)  define  an  operator  A:  H\C)  - — > Z and,  for 


fixed  2,  a subset  of  X,  supposed  non  empty: 

I,  . X j j 


(3.2) 


The  null  spaces  of  the 


operators  A and,  D are  given  by; 

I <r  Kc,*>h^  < 

to (Dij)  --  ^ xe  //Vo  / ojsjxe^fc)  I x--<W./ 


Their  dimensions  are  equal  to  (2n)  and  to 


one,  respectively,  so  that 


> 1. 


(3.3) 


A/  C ^ ) r)  N [f\y  ^/provided  n 
Hence,  J , the  minimum  problem: 

(b  l x'V(z)J&fz 

tel;  JC 

I*J  xte),  tCj  < to:'/X>0, 

~ x Hi  ± n l 

? * (*>/.-,*«.  Sr/,  —.,1:4  .)  5 ft"' 
h“  * “,“’ue  SOlUtion  ^ O C #Vcj  . for  *,  as  long  „ n>  ,, 

This  unique  element.^  f i,  herein  -a, led  the  He  mite  interpolating 
closed  spline  corresponding  to  f H'(c),  fly  3. 

The  n)-  dimensional  subspace  ^htc  //  C^“  J 
closed  splines  is  defined  by: 


of  the  Hermite 


uc 


C ft 

*1  i fi  If  tcJ  I * fy  * > e=  ° ; -bx e A/(7?> 1„  j 


(3.4) 


any  ^ , there  exist  a unique  element  G"£  S such 

that  ft<T-Z  . It  has  the  extremal  properties  described  by  equati 
(2.fl)  and  (2.9)  /T-2J. 


ons 
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Existence  and  uniqueness  of  the  solution  of  the  minimum  problem 
(3.o)  (.and,  hence;  of  the  problem  is  thus  proved. 


i'roofs  of  the  characterization  and  properties  of  the  Hermite  closed 
interpolating  splines  discussed  in  the  previous  paragraph  hinge  on 
the  following  theorem  b 7: 

Theorem  1 
‘ 1 T 

^“6  closed  Hermite  interpolating  function  corresponding 

to  1 , Kc,  Ki  > t-  t fc;  ) if  there  exist  (2n)  coefficients  and  A' 

such  that: 


Me  t / \i  Ki  )e  j a^}J 


1 


(3.5) 


Here  D.  denotes  the  adjoint  of  Dfl  and  L J 
1 7 


± 


// 


the  orthogonal  subspace . 

The  demonstration,  with  few  variations,  is  substantially  the  same  as  that 
employed  in  [ l J and  only  the  essential  points  will  be  detailed. 

The  definition  of  adjoint,  eq.  (3.5)  with 

$ ~/  1 ‘ 

A / = C -IJ 

J 


A;-C-0U:  ■ A :=c-0r  A; 


and  equation  (3.1)  lead  to: 


/ 


U.<r,  b,  x > --  k.  (-if  [ k x (t;J -A[-  X '[<■)] 

H 4 = / 

V x 6 14  Yc  J 


from  wnict.  one  deduces  that.OJ" 

a)  the  conditions  imposed  by  theorem  1 reduce  to  the  only  requirement 


«...  £/d,-KVv|7ie>  //*(  orthogonal  t 


o unity  and  this  leas  to: 


v / 

.>  A ■ - o 

i ^ 1 


(3.6) 
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b)  Characterizing  S'  amounts  to  finding  a function  fe  H“(C.) 
which  verifies  the  identity 

i (-')’>  L^i  t'fcJl  - *-f,  D,  *>y0 

VteH'Cc) 

subject  to  the  conditions  (3.6)  for  then  one  has: 


(3.7) 


• >-1/ 
with  Y an  arbitrary  constant  and  such  an  expression  for  (5*  wi] 

°o 

contain  (2n)  arbitrary  independent  parameters. 

Mac-Laurin  developments  of  X ) and  X ' with  the  rest 
expressed  in  integral  form  yields: 

X (T.  j-  £ *(<%  + (b  ‘x%)dz 

r~c  /!  Jc  {^-0! 


y zb‘lh(o)  , ^ Cr; -i)T 


xbih't-R-LL 

h 0-oj 

Account  for  the  identity 


(3 


tr-Zcb‘ c-')*  (r  - X + 


apply  the  binomial  formula  to  J 


(?.•  - *) 


and  define : 


Yl  - / 


L 1 yri  j ( m-/  JJ 


-»!>(? 


(3.8) 


to  obtain:  £ since  0)o  = 0 upon  eq.  (3.6)  and  6/0  = 0 by  definition  7: 


< t.  *>  «=  £ [ l *fc)  * £ 


H 


/ 

= £.  [d; +*/,■  ] {"'to)  + 


i - / 


y j 


(3.9) 


Ij?  [ p™?'.  jinxt, 

r hz‘  Q-  i)J  r<* - z)i  * f 


(3  -2  )J 


with: 


v / $-/  a - 2 ' — • 

(?/r):E  yl»  fe  . ]>Z  (-if  el 

0 ' /C.  . I I r>  J ■?- 


z' 


'■*  jv 


(3.10) 


H 


— 9-1  8-2 

t((r)r  2>_  i,-  (r.--zj_  . JTJ 


t 


<J  - z)  J 4 


- C 


(->)*  J I hzL 

<j! 


Only  the  first  term  on  the  right  hand  side  of  equation  (3.9)  must  be  further 
transformed  to  put  it  in  the  form  of  a scalar  product  in 

'*  jk  ~~  2 ' : - » •**  proven  i2i  (1  J t tne  (absolute)  con- 


tinuity  of  * frj  for  implies  that 


7 


-L  / A:  t 7;  / r a,  j r <r  tr  r\ 


.»  / / b y 


' I U-  ) 


, r,  U°(c  1 


f 7 VI  ) 
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where : 


t:' 

if  (?)  = Z-  Y: 

, SI  0 


if 
0 * 


with  the  coefficient  ^ such  that 

0 

z 


«/ 


Cj+t  t/  _ 

= C-0 


I 6 t 6 (Sj~/ 

Substitution  of  eq.  (3.11)  into  eq.  (3.9)  shows  that: 


(3.12) 


(3.13) 


Hence,  on  account  of  eq . (3.7),  the  following  characterization  of 
is  obtained: 


± (±102  ]..k  B . z± 

(j  -i)J  J j'0  r’tt  jl 


•:  -.re,  upon  eqs . (3.10)  and  (3.12): 


|^+J  ~Yj  ^ (3-15) 

o±}±yi 

• I 

Of  the  (2n  + q)  coefticients  ( B - , ^ l ) appearing  in  the 

/a  \ ' *j  r y ' 

expression  (3.lM  for  5T  onIy  (Pn)  are  independent. 

They  are  indeed  related  by  eq.  (3.6)  and  by  the  last  (q-l)  equations 
(3.15)  which,  on  substituting  ft's  (for  J O ) from  eq.  (3.13) 
can  be  rewritten  as:  ” 


2—  I Bttj  - t +, 

>=0 L ThOJ 


C-o 


Cj-r'tTj 


Lo 

t / ’ — 1 


(3.16) 


/ ir  t ± 6j 


These  equations  for  ( 3 5:  ^ ^ a - / ) express  the  vanishing  of  the 

(tt  I 1 

discontinuities  of  (5~  and  of  its  first  (q-3)  derivatives  at  the 

closure  point,  as  the  following  Lemma  shows: 


Lemma 


The  function  ( Z ) and  its  fi  rf;t  ( q - 3 ) derivaties  are 


continuous  on  C " 


Proof 

Given  the  properties  of  polynomials  and  of  the  functions  fe" - Z"i)  the 

statement  of  the  lemma  needs  to  be  proved  only  at  the  closure  point  . 

In  the  interval  J Zt  = Z [ eq.  (3.1U)  reduces  to 


In  the  interval  J 2" £ 

rl'h:  ^-6  - 


'cz)=  ?-c  P . zL  , L. ill , M 


i f-  Z- 
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In  the  interval  3 t[  1 1 <■  4 (Z~  for  any  ( i ) so  that, 

on  account  of  eqs . (3.10),  eq.  ( 3 . 1 U ) yields: 


V' 


I'2 


o 


OIjLZ—  - 
( 


a r* 

The  discontinuities  of  the  (q  + k)-th  derivative  at  the  closure  point 

are  therefore  given  by 

- *L°"k,(o')-^kJ(ry- 

= - -51  / JWtit' . t C-i)  f > I 

L-  , - 


/ / t / / ■' 


or,  with  k + 1 * q - r: 


o<k 


(2^-t  -<) 


“ PR 


X.  6* 


/CM 

^ u) 


- o 


L 


11 


(hOJ 


2it£  sj- / 


*)**■*<},-, 

+ ( - 0 

j! 

(3.17) 


and  comparison  with  eq.  (3.L6)  comnle  + es  t.he  nn-of  of  rbr>  i^pina. 

The  hermite  closed  spline  functions  t>  ( 'Z  ) corresponding  to  the 
set  j are  obtained  by  integrating  equation  (3.U)  7 - times  and  using 
Krr.i  trary  constants  to  endorse  the  continuity  of  s (Z  ) and  of 
its  first  (q  - 1)  derivatives  at  the  closure  point,  sc  as  to  make  A <c  H'&-  ) 


Vli'  1 u-:-  ,;j  characterization  of  s (Z)  given  by  eq.  (2.2), 

’ " ' ' y h :ks  ■ ha  fcJ  - - • ut.  t>  oi 


On 
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(2q  - r - l)-the  derivatives  at  the  closure  point  are  still 

given  by  eqs.  (3.17)  with  >z  now  ranging  from  2 to  (2q  - l). 

Hence  imposing  the  continuity  of  s (T)  and  of  its  first  (q  - l) 
derivatives  amounts  to  the  statement  that  eqs.  (3.16)  must  hold  also 
for  (r)  ranging  from  the  value  (q)  [ corresponding  to  the  condition 
'b  ^ the  value  (2q  - l)  ^"corresponding  to  the 

condition  £ ^(Pi)-t?J.  Equations  (2.U)  are  thus  proved. 


This  completes  the  proof  of  the  statements  made  in  paragraph  (2). 


In  particular,  equation  (2.3) 
(2.7)  is  the  definition  (3-M 


is  the 
of  the 


condition  (3.6)  and  equation 
subspace  S//c  • 
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